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ABSTRACT 



Context. Scattered light images of circumstellar disks play an important role in characterizing the planet forming environments around 
young stars. The characteristic size of the scattering dust grains can be estimated from the observed brightness asymmetry between 
the front and back side of the disk, for example using standard Mie theory. However such models often overpredict their brightness 
by one or two orders of magnitude, and have difficulty explaining very red disk colors. 

Aims. We aim to develop a dust model that explains simultaneously the observed disk surface brightness, colors and asymmetry in 
scattered light, focussing on constraining grain sizes. 

Methods. We use the 2D radiative transfer code MCMax with anisotropic scattering to explore the effects of grain size on synthetic 
scattered light images of circumstellar disks. We compare the results with surface brightness profiles of the protoplanetary disk HD 
100546 in scattered light at wavelengths from 0.4 to 2.2 micron. 

Results. We find that extreme forward scattering by micron sized particles lowers the effective dust albedo and creates a faint and red 
disk that appears only slightly forward scattering. For the outer (> 100 AU) disk of HD 100546 we derive a minimum grain size of 2.5 
micron, likely present in the form of aggregates. Intermediate sized grains are too bright, whereas smaller grains are faint and scatter 
more isotropically, but also produce disk colors that are too blue. 

Conclusions. Observed surface brightness asymmetries alone are not sufficient to constrain the grain size in circumstellar disks. 
Additional information, such as the brightness and colors of the disk are needed to provide additional constraints. 
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1. Introduction 

Protoplanetary disks are thought to be the main sites of planet 
formation, and spatially resolved scattered light images play 
an important role in revealing what happens before and dur- 
ing the birth of planetary systems. With superior diffraction- 
limited resolution compared to thermal emission at longer wave- 
lengths, scattered light images in general - and in particular 
those from the Hubble Space Telescope - have provided us with 
some of the most detad ed images of protoplanetary disks so far 
(e.g. Grady et al. 2001). Recent advances in adaptive optics and 
data reduction techniqu es have greatly improved the efficacy of 
ground-based data (e.g. Thalmann et alJl20lol : [Hashimot o et alJ 
1201 U lOuanz et al.l 1201 lb . while future observatories such as 
ALMA will achieve similar spatial resolution at (sub)millimeter 
wavelengths. 

Scattered light images offer important information on the 
overall geometry of protoplanetary disks and contain possi- 
ble signposts of planets. Large scale spiral structures have 
been observed, some trigg e red by a close encou nter from out- 
side dClampin et all 120031: lOuillen et ail 120051) . others more 
likely by a (planetary) perturber from within (e.g. iGradv et"al] 
120011) . Annular gaps and inner holes are also observed, and 
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are likely carved by orbiting protoplanets (IGradv et alJ [2007; 
iThalmann et"al]|2010l:lHa"shimoto et al.ll201 lb . 

Apart from the macroscopic signposts at the end of the (gi- 
ant) planet formation process, scattered light imaging can also 
shed light on the earliest stages of planet formation: the growth 
of dust grains. The main diagnostics for this are: 

1. Brightness asymmetries between near and far side of the 
disk. In an inclined disk, the near and far side of the disk 
surface are seen under a different angle with respect to the 
star. If the dust scatteres anisotropically, a brightness asym- 
metry is observed. Larger dust grains are more forward scat- 
tering, and will show a stronger brightness asymmetry. The 
observed values for the asym metry parameter g in differ- 
ent disks are betw een -0.15 dArdila et alJl2007l) and -0.8 
dPinte et al.ll200 8b). ranging from nearly isotropic to slightly 
forward scattering. 

2. Disk surface brightness or dust albedo. The fraction of 
scattered light at a specific wavelength. In general, disks 
are observed to be much fainter that predicted on basis 
of their grain size. A sample of Herbig stars imaged by 
Fukagawa et al.ld2010l) shows fractional luminosities in scat- 
tered light on the order of a few percent or less. This is 
much too low for non-Rayleigh scattering particles, which 
typically have an albedo of ~ 0.5. In debris disks, which 
are optically thin, albedos on the order of 0.05. ..0.1 are de- 
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Fig. 1. Visualization of anisotropic scattering properties of different dust species. Displayed are: a small particle in the Rayleigh 
limit (2na < A) in blue, an intermediate sized particle (2na ~ A) in gray and a large particle (2na > A) in red. The angle between 
the disk surface and incident stellar light is 3 degrees, equal to the flaring angle minus aspect ratio of the surface at 200 AU. The 
phase functions at 0.6 pm are plotted in polar coordinates, and normalized such that surfaces are equal (square root of intensity). 
The arrows indicate the direction and magnitude of light scattered towards the observer. 



rived from observations dKalas et alj|2005b iKrist et aDl20 1 Ot 
iGolimowski et al.ll2.Ql lh . while an albedo of 0.5 is predicted 
on basis of the grain radius of ~ 1 pm derived from asymme- 
tries in the disk images. 
3. Disk color, caused by the wavelength dependence of the 
albedo. Dust grains in the interstellar medium have a grey 
scattering colo r at optical to near-infrared wavelengths (e.g. 
IWhitnevlll995l), whereas more evolved objects such as de - 
bris disk s (e.g. IGolimowski et all 12006: De bes et alJ l2008). 
comets (Jewitt & Meech] 119861) and Kuiper Belt objects 
(e.g. iLuu & Jewitt 19961) can have redder colo rs, indicat- 
ing th at larger than ISM grains are present (e.g. IWolf et al] 
2012h. Most protopl anetary disks have grey colors (e.g. 
Fuka gawa et alJ201 |). though some di sks exhibit redder col- 



ors: H P 141569 dClampin et alJl2003h. G G Tau dKrist et alJ 
l200l and HD 100546 (Urdila et al]l2007l) . 

To extract the grain size from scattered light images, these 
three aspects (moderate brightness asymmetries, low albedos, 
neutral to red colors) have to be modelled simultaneously. This 
is challenging, because small dust grains have low albedos but 
very blue colors. Large grains have neutral to red colors but also 
high albedos. We call this "the color and brightness problem of 
scattered light images". 

To solve this problem, we have studied the effects of 
anisotropic scattering on the observed brightness and colors. 
Particles that are large enough become extremely forward scat- 
tering. The bulk of the scattered light is concentrated in the for- 
ward peak - which is outside the range of observed angles for a 
disk that is not observed edge-on (Fig.[TJ. The resulting few per- 
cent are scat tered towards the observer , resulting in a low 'effec- 
tive' albedo dDullemond & Nattal2 003). Because the observer of 
an inclined disk sees only a part of the phase function outside of 



the forward scattering peak, the disk image does not appear to 
have a strong br ightness asymme try, and can even appear back- 
ward scattering dMin et al.ll2010h . We will describe this effect in 
more detail in section [2] 

We have applied this strong forward-scattering dust model 
to the well-studied protoplanetary disk of HD 100546, for 
which scattered light i mages are available from 0.4 to 1.6 pm 
lAugereau et alJ d2001l) : lArdila et all d2007), and we present a 
new image at 2.22 pm in section [3~TI After constructing a ge- 
ometric model of the disk in section 13.31 we derived the grain 
size from the observations in sections [3~4l and 1331 We will also 
discuss how the derived grain size compares to those of other 
size indicators in the discussion. 



2. Scattering by protoplanetary dust 

In this section we will briefly discuss the basis of the scatter- 
ing p r operties of protoplanetary dust dMid [l908: Va n der Hulsj 
[T957l lBohren& Huffman! [1983) and its influence on the ob- 
served brightness and color of scattered light images. 



2.1. Single scattering albedo 

The fraction of light that is scattered by a dust grain in all di- 
rections is given by the single scattering albedo, 10. For particles 
much smaller than the wavelength the albedo is very low since 
most light is absorbed by the particle. When the size of the par- 
ticle increases, the single scattering albedo also increases. For 
particles much larger than the wavelength of incident radiation, 
the total extinction cross section of the particle is comprised of 
three components: absorption (C a b s ), reflection/refraction (C re f), 
and diffraction (Cdiff). 
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Fig. 2. Real versus effective albedo of particles of different sizes. The effective albedo (solid line) is defined as the fraction of light 
scattered outside of the forward ten degrees times the real albedo, see text (dotted line). The green area denotes the range of observed 
wavelengths by the Hubble scattered light images. The small, medium and large panels correspond to a particle size of 0.08 //m, 
0.25 yum and 2.5 jum. 



The absorption and reflection/refraction cross sections come 
directly from geometrical optics and sum up to the geometrical 
shadow of the particle. The diffraction cross section comes from 
the distortion of the wavefront caused by the particle and is also 
equal to the geometrical shadow of the particle. It is than directly 
derived that the single scattering albedo of large particles is, 

C sca t C re f + Cdiff > Q 5 ( 1 ) 



C e 



Cabs + C re f + Cdiff 



2.2. The effective albedo and the color of scattered light 

However, the single scattering albedo is an angle-integrated 
quantity, and care should be taken when scattering becomes in- 
creasingly more anisotropic. In the limit of Rayleigh scattering - 
where particles are much smaller than the observed wavelength 
- the phase function (Fn(8)) of scattering is close to isotropic 
(Fig [3]). For particles larger than the wavelength the scattering 
becomes anisotropic and heavily peaked towards the forward di- 
rection when the particle size increases. 

Light which is scattered in the forward direction by a parti- 
cle in the surface layer of the disk is scattered into the disk, and 
will not be detected in scattered light images (see Fig. [1] large 
particle). Therefore, it is clear that for particles that scatter pre- 
dominantly in the forward direction the effective scattering cross 
section, i.e. the amount of light that is actually scattered into our 
line of sight, is much smaller than would be expected from the 
angle integrated single scattering albedo. Therefore we define 
an effective albedo which only takes into account the part of the 
phase function that can actually be observed between angles 8\ 
and 02'. 
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Fig. 3. Phase function of a large dust grain (2.5 fim, red line), an 
intermediate size dust grain (0.25 /mi, grey line) and a small dust 
grain (0.08 //m, blue line) at 0.6 //m. The green area marks the 
range of observed angles for HD 100546, between 34° and 126°. 
8 = 0" refers to forward scattering and 8 = 180° to backward 
scattering. 

For large particles this effective albedo can be smaller than 0.5 
when the phase function is sufficiently peaked towards the for- 
ward scattering direction. 
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Fig. 4. PSF subtracted NICMOS images at 1.6 and 2.22 /im. The black region in the F160W image is masked out due to a strong 
diffraction spike. 



In Fig. |2]we plot for three different grain sizes the effective 
albedo for arbitrareljQ G\ = 10°, G 2 = 180° together with the total 
albedo. For the small grains the scattering is quite isotropic (see 
Fig. [3J and the effective and total albedo are almost the same. 
The color of scattered light for these grains is always grey to 
slightly blue. 

The phase function of larger particles is more anisotropic, 
making the effective albedo quite low, even below the limit for 
the single scattering albedo (Eq.[T]l. When going to longer wave- 
lengths the phase function becomes more isotropic and the ef- 
fective albedo increases. This increase of effective albedo with 
wavelength will be reflected in the observed color of the scat- 
tered light images, producing reddish colors. 

3. Application to HD 100546 

To test if exteme forward scattering can explain the observed 
faintness of protoplan etary disks relative to their host stars 
dFukagawa et al.l uOlOh . as well as their color indices and ob- 
served brightness asymmetries, we have compared our dust 
model to scattered light images of a well-studied Herbig star, HD 
100546. Its disk has been imaged by the Hubble Space Telescope 
in sc attered light over a broad wavelength ran ge of 0.4 /zm to 1 .6 
fim dAugereau et alJl200lt lArdila et all 120071) . We will extend 
this wavelength range using a new NICMOS image at 2.2 /zm. 
This large wavelength range allows us to compare both bright- 
ness and colors. We will describe these observations in the next 
section. 

3.1. Observations 

We took HST/NICMOS coronagraphic images (0.3" radius im- 
age plane obscuration, camera 2 image scale 75.8 mas/pixel) of 

1 Cutting out the forward 10 degrees is sufficient to illustrate the main 
effect of not observing the forward scattering peak 



HD 100546 and the PSF reference star, HD 109200 with the 
F222M filters (central A = 2.22 /zm, A/AA ~ 10) on 16 March 
2005 as part of Program GO 10167 (PI: Weinberger). The ob- 
servations include long exposures with the stars underneath the 
coronagraphic spot for high contrast imaging at two different 
spacecraft orientations and direct images of both stars outside 
the coronagraphic hole with short exposures for point source 
photometry. The instrumentally calibrated and reduced images 
discussed in this paper were created from the raw NICMOS 
multiaccum exp osures following the pr ocessing methodolgy de- 
scribed by §3 of lSchneider et al.l ([2005) and references therein. 

For photometric analysis, each calibrated direct image was 
used separately to independently determine the total photometry 
of the star and empirically determine the uncertainties in each 
filter band. The three images for each star and in each filter were 
located at different positions on the detector. 

We used a median combination of the three dither points to 
create a final image of each star to derive a ratio for scaling and 
for photometry of HD 100546. We used a 20 pixel radius cir- 
cular aperture to determine the photometry. The background in 
the images is zero, so no background annulus was used. The in- 
dividual dither points were used to get a rough estimate of the 
uncertainty in the ratios and photometry. At F222M we measure 
a total F v =5.4+0.1 Jy, and a scaling with the PSF reference of 
0.89+0.04. 

In order to determine the best subtraction we minimized a 
chi-squared metric on a region of the target image dominated 
by the star's diffraction spikes. We assumed that good subtrac- 
tion of the diffraction spikes in a region uncontaminated by the 
disk corresponde d to the best subtractio n of the PSF within the 
region of interest (Schn eider et al.ll2001h . We iteratively created 
subtractions for combinations of scaling and pixel offsets until 
we found an image that produced the lowest chi-squared mea- 
sure. We searched within ±1 pixel to find the best x and y pixel 
offsets. 
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To quantify the systematic effects on the photometry, we re- 
peated the subtractions varying the PSF scalings and offsets by 
+ 1 cr from the minimum chi-square solution found above. Using 
an elliptical photometric aperture matched to the estimated in- 
clination of the disk and that extended from between (X'5 to 4", 
we found the standard deviation in the disk flux densities from 
this suite of subtractions. In F160W and F222M, the total mea- 
sured flux of the disk at both spacecraft orientations matched 
to within the uncertainties. We then propagated this systematic 
uncertainty into the total uncertainty in the flux density of the 
disk per pixel. Subtracted images were then geometrically cor- 
rected for the slight optical distortion of the NICMOS camera 2 
at the coronagraphic focus. We used the x-direction pixel scale of 
(X'07595/pixel and the y-direction pixel scale of O'.'07542/pixel 
to create an image with pixels that have the y-direction plate 
scale in both directions. The geometrically corrected images 
were rotated about the position of the occulted star to a com- 
mon celestial orientation using the rotation centers given by the 
flight software in the raw data file headers. The final result is 
show in figure |4] 

Additional observations of HD 100546 with NICMOS were 
perfo rmed as a part of the GO program 9295 (Augereau et alJ 
1200 lb . Observations with NICMOS in the F160W filter were re- 
covered from the archive, reduced in the same manner as the 
medium bandwidth image, and median combined. Archival PSF 
reference stars for the images were subtracted from the target 
observations. The F 160W reference was HD 106797 as used in 
Auge reau et al.1 (l200lh but we followed the above procedure for 
subtraction as with the F222M filter. Figure|4]shows the resulting 
PSF subtracted NICMOS images of the HD 100546 disk used in 
this work. 

HD 100546 was also observed with the HRC coronagra ph on 
ACS in the F435W and F814W filters dArdila et alj|2007b . The 
final reduced images were kindly provided to us by D. Ardila. 
For the ACS observations, no independent estimate on the un- 
certainties in flux scaling for the disks was performed. We as- 
sumed that for both instruments this affect was ~5%, comparable 
to what we calculated for our medium bandwidth filter observa- 
tions. 



3.2. Surface brightness profile and disk albedo 

The surface brightness profiles at 0.4, 0.8, 1.6 and 2,2 fim are 
shown in figure [5] They were constructed by taking the median 
surface brightness over an elliptical annulus whose shape cor- 
respondspl to a disk at a position angle of 145° and inclin ed at 
46° - the average in clination inferred by Pa ntin et all (l2000l) and 
Ardil a et all (|2007). To be able to compare the surface bright- 
ness profiles of images taken at different spatial resolution, the 
width of the annuli is taken to be 0.15", the diffraction limit of 
the longest wavelength filter. For the uncertainties we used local 
estimates at the major and minor axes of the standard deviation 
of counts within smaller apertures. Missing data were treated as 
NAN values and were not included in the median. 

To measure the reflectivity of the disk observed in each filter, 
we define a geometric albedo Wdisk for the disk, which is the 
fraction of (star)light scattered towards the observer at a specific 



2 Although the flux coming from a certain distance from the star devi- 
ates from an ellipse if the disk is flared, we don't expect this to influence 
our analysis as long as we compare observed and model images in the 
same way. 
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Fig. 6. Geometric albedo Wdisk of HD 100546 in the different 
filters. The geometric albedo is a measure of the fraction of 
light scattered by the disk (see eq. and text for definition). 
Displayed are both observations (diamonds with systematic and 
measurement errors) and the same models as figure |5j A large, 
intermediate and small grain (red, gray and blue solid line re- 
spectively). A model with intermediate sized grains with an 
albedo close to one is also displayed (gray dotted line). 



radius r. It is defined as: 

SB(r) 

where SB(r) is the surface brightness of the disk in scattered 
light, F„ is the flux density from the star and F mnf , r is the flux 
density from the inner (<1 AU) disk, both measured from the 
unresolved Note that this 'disk albedo' Wdisk depends on 

both the dust albedo and the geometry of the disk (i.e., its flaring 
angle). 

To be able to quickly compare disk albedos at different wave- 
lengths for both the data and models, we define a geometric 
albedo per filter a)^ at a typical radius of 2". We obtain com 
by fitting a powerlaw of the form a»ait ■ ( r )~ q to a»disk (f) between 
1.0"and 3.5". When plotting this disk albedo versus the wave- 
length of the filter (Figure|6), we obtain a diagnostic tool that de- 
scribes the disks geometric albedo and colors. Throughout this 
paper, we will apply this same procedure to the synthetic model 
images as well. 

As can be seen from figure|6] the disk is relatively faint com- 
pared to the star and quite red over the entire wavelength range. 
Even though HD 100546 is one of the brightest disk to be ob- 
served in terms of absolute flux, it reflects only about a percent 
of the incoming light in the near infrared, and about a factor 
of 4 less in the bluest filter. Any dust model invoked to explain 
scattered light images must be able to explain both the relative 
faintness and red colors. However, as mentioned before, the ge- 
ometric albedo is a product of both dust properties and disk ge- 
ometry. Therefore a geometric model of the disk is necessary to 
isolate the effect of the dust albedo. 



3 Note that at near-infrared wavelengths, the thermal emission from 
the inner disk is much brighter than that of the star. The light that we 
observe in the longest filters is therefore dominated by scattered inner 
disk light, rather than scattered starlight. 
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Fig. 5. Radial surface brightness profiles of scattered light images of HD 100546 (diamonds, grey area indicates systematic and 
measurement errors) and disk models with: A large dust grain (2.5 yum, red line), an intermediate size dust grain (0.25 /vm, gray line) 
and a small dust grain (0.075 yum, blue line). 



3.3. Disk model 

To interpret the scattered light images of HD 100546, we will 
use the 2D radiative transfer coddj MCMax (iMin et all [2009) 
with anisotropic scattering for the dust as described in section|2] 
This code has been succesfully applied for modelling sever al ob- 
servables of protoplanetary dis ks (e.g. Verhoeff et al.ll20l"il) and 
includes anisotropic scattering (IMin et al 120121) . 

In constructing the disk model, we do not have to start from 
scratch, as the disk geometry for HD 100546 is already well con- 
strained by previous observations and mod els. We will use the 
SED fit presented by iMulders et al.l (1201 ll) as a starting point 
for our model. The model features a depleted inner disk start- 
ing at 0.25 AU (iBenistv et al.ll2010tlTatulli et al.ll201lh which is 
fairly small (<1 AU. iPanic et al.ll2012l) . has a large empty gap 
( Bouwman et al] 120031) and a massive outer disk star ting from 
13 AU to 400 AU dGradv etaDl2007llPanic et al.ll2010h . 



The vertical density structure of the disk is described by 



), where the 



fitted scale height 

1.3 



a Gaussian (p c 

is parametrized as H v = 0.04AU(r/AU) 1J , consistent with 
Mul ders et al.l (1201 ll) . The reason to deviate from solving the 
vertical hydrostatic structure is that we want to study the effect 
of the dust properties in a fixed geometry, whereas in a hydro- 
static disk the geometry changes with dust properties through the 
temperature. For this purpose it is therefore easier to use the pre- 



4 Since both scattering and absorption take place in 3D, this model 
can also be referred to as an axisymmetric 3D radative transfer code. 



Parameter 


Value 


M disk 


0.0001 M Q 


P 


1.0 


inclination [°] 


46 


Rin [AU] 


0.25 


Rgap,i„ [AU] 


0.3 


Rgap.ou. [AU] 


13 


Rou, [AU] 


350 


f. 

Mnner 


0.05 


HiAU 


0.04 




1.3 



Table 1. Disk parameters for the geometrical model. The surface 
density profile is defined as S(r) oc r~ p and scaled to the total 
disk mass. The scaleheight is defined as H p {r) = H\au^- The 
gap ranges from R gap>in to Rg ap , out 



scribed vertical structure based on the hydrostatic model, which 
provides an equally good fit to the SED (Fig[8]). 

The inner disk (< 1 AU) is not directly probed by our images, 
but its thermal emission is a factor of a few brighter than the star 
in the near-infrared. Light scattered of the outer disk surface in 
the longest filters is therefore domi nated by scattered inner disk 
light, rather than scattered st arlight (Pinte et al]|2008al) . We take 



the dust composition from Muld ers et al.l (1201 ll) for the inner 



disk (table 0, and keep this composition fixed if we vary the 
composition of the outer disk, such that the light illuminating 
the outer disk remains constant. A first confrontation with the 
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1.0 



53 0.5 



0.6 



0.4 



0.2 



0.0 
0.01 



Grain 


size \pm ] 


composition 


large 


2.5 


Solar 1 


intermediate 


0.25 


Solar 1 " 


small 


0.08 


Solar 1 


inner disk and wall (<20 AU) 


0.4 


ISM 1 



0.10 



.00 
a 



10.00 



Table 2. Dust properties of our best fit model. Opacities are 
calculated assuming a grain size distribution f(a) oc cr 3 - 5 
with typical size a with a width of 0.25 dex. f Composition 
of ISM: 13.8% MgFeSiC-4, 42.9% MgSi0 3 , 38.3% Mg 2 Si0 4 , 
1.8% NaAlSi 2 6 , 5% C, Solar: 12 % MgFeSi0 4 , 12 % 
MgFeSi 2 6 , 12 % Mg 2 Si0 4 , 12 % MgS i 7 O fi , 15 % FeS 40 % 
C. Optical constants are from : Silicates dDorschner et al.l fl995; 
iHenning & Stognienkol Il996t iMutschke et al.1 1 19981) . Carbon 
1 00.00 (Preibi sch et al.lll993l) . Troilite (iBegemann et al.l ll994). 



Fig. 7. Asymmetry parameter g as function of particle size at 
0.6 pm, calculated in two different ways. The solid line uses 
the formal definition g = (cos 0), the dotted line is a Henyey- 
Greenstein function fitted to the phase function in the observed 
range of scattering angles. The observed range of values of the 
asymetry parameter is indicated by the gray area. 



observed scattered light images shows that using the same com- 
position for the outer disk does not work well. It overpredicts the 
scattered light flux by a factor of 10-30, see figure[6]gray dotted 
line. 

However, the albedo of this dust is very high, close to one. 
We first explore the effects of lowering the albedo through the 
dust composition to reach a lower surface brightness. We switch 
to a different composition which has a particularly low albedo 
of ~0.5. This is is also the theoretical limit in the geometri- 
cal optics regime (see section |2| and in practice, a lower real 
albedo can only be achived by changing particle size (See Fig 
12), which we will do in the next section. For an albedo of 0.5 we 
use a dust co mposition derived from a condensation sequence, 
as described in lMin et~ai1 (1201 ll) . This lower albedo reduces the 
amount of scattered light significantly, but still overpredicts the 
observations by an order of magnitude (Pig. O dotted line). The 
final parameter we have to tune to obtain the low observed sur- 
face brightness is the grain size, which we will do in the next to 
sections. 

3.4. Brightness asymmetry 

To constrain the characteristic grain size in HD 100546 we look 
first at the observed brightness asymmetry of the disk. As already 
explained in section |21 particles smaller than the wavelength of 
light (27Tfl v < A) scatter isotropically, whereas large particles be- 
come increasingly more anisotropic. This anisotropy can be de- 
scribed by the asymmetry parameter g, which is the expectation 
value of the cosine of the scattering angle (g = (cos8)), a param- 
eter between and 1 that increases with grain size (Fig [7). 

In reality we can only sample a limited part of the phase 
function (34 to 126 degrees in HD 100546, for a 46°inclination 
and a 10° opening angle), and g can not be calculated directly. 
The phase function is therefore oft en approximated by a Henyey - 
Greenstein function /hg(,?hg,#) dHenvev & Greenstein 1941), 
which is parametrized in such a way that the expectation value of 
the cosine of the scattering angle is equal to guG- This way one 
can derive the asymmetry parameter with only partial knowledge 
of the phase function. 



However, care should be taken by deriving grain sizes in 
this way. When grains become extremely forward scattering 
(g> 0.5), the forward peak of the of the phase function caused 
by diffrac tion moves out o f the observed range of angles (See 
section El MinetaU 12010). The observed range caused by re- 
flection/refraction is not necessarily forward peaked (Fig. [3), and 
the derived value of the asymmetry parameter geG can become 
much smaller than the real asymmetry parameter g (see figure 
®. 

The disk of HD 100546 shows only a small brightness 
asymmetry of j?hg = 015 at .4. ..0.8 pm, with maxima up to 
guc — 0.23 (lArdila et al.ll2007h . Comparing this number to the 
real brightness assymetry g yields a grain size of a ~ 0.08 pm, 
similar to a = 0.1 /mi derived by Pan tin et al.1 yOOO) based on 
gHG - 0.2. However, taking into account the limited part of the 
phase function observed yields a second set of solutions with 
grain sizes larger than one micron which come close to the ob- 
served asymmetry, see figure [7] 

3.5. Disk color and brightness 

Although small particles (a = 0.08//m) provide a good fit to the 
observed brightness asymmetries, they can not explain the disk 
color. They scatter in the Rayleigh limit (27ra < A), and hence the 
scattering efficiency drops off as /T 4 . Such a decreasing albedo 
with wavelength (Fig.|2) will give rise to very blue disk colors, 
and it becomes clear that these small particles might explain a 
low albedo at one wavelength, but never the observed disk col- 
ors (Fig. [6] blue line). As discussed before, intermediate sized 
particles (2na ~ A) have a high albedo > 0.5 and overpredict the 
disk brightness (Fig. [6] gray line) and also appear too forward 
scattering (Fig|7). 

Larger particles (2na v > A), on the other hand, do a much 
better job in explaining the disk colors and brightness. As de- 
scribed in section|2]and shown in figure|21 forward scattering de- 
creases the effective albedo towards shorter wavelengths, though 
the real albedo stays above the theoretical minimum of 0.5 (Eq 
[TJ- A particle of 2.5 pm fits both the colors and brightness of the 
disk (Fig. [5] and [6] red lines), though its brightness asymmetry is 
a little too high. We will discuss this further in section l4~3l A par- 
ticle with a different composition that has a higher real albedo 
would need to be larger to reduce its effective albedo to the same 
level. For example, a particle of different composition that has a 
real albedo of 1.0 needs to have a size of 10 pm to explain the 
observed brightness. 
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X [yiim] X [l J - rn ] 

Fig. 8. Observed SED of HD 100546 plus photosphere (grey lines). The solid lines are two model SEDs with our best-fit aggregate in 
the entire outer disk (left) and with the aggregates only outside of 20 AU (right). The dotted and dashed lines denote the contributions 
to the SED of scattered light and thermal emission, respectively . 



4. Discussion 

We have shown that particles larger than 2.5 micron can explain 
the red colors and low brightness in the disk of HD 100546. 
The question arises whether such large grains are consistent with 
other grain size indicators in the SED, and why we overpredict 
the brightness asymmetry, which we will discuss in the next sec- 
tions. 

4.1. The 10 micron silicate feature 

One grain size indicator in the SED is the 10 yum silicate feature, 
which indicates the presence of micron-sized grains or smaller. 
However, at a wavelength of 10 /im particles larger than 2.5 yum 
scatter extremely efficiently because they are in the resonance 
regime (27ra ~ lOyum). In addition, for particles this large com- 
pared to the wavelength of 10 yum, the absorption and scattering 
coefficients add up to the geometrical shadow of the particle, as 
is the case for particles in the limit of geometrical optics (Section 
0. Hence the silicate feature is imp rinted inversely on the scat- 
tering efficiency at this wavelength ( M in et al. 2004). The com- 
bination of both effects supresses the silicate feature (Fig. [8]). 

However, it should be n oted that the sil icate feature arises 
mainly in the disk wall dBouwman et alj2003h . whereas the scat- 
tered light images trace a region outwards of 0.5"(~50 AU at 
a distance of 103 pc). Radial variations in the dust properties 
have been observed in scattered light (e.g. lOuanz et al.ll201ll 
Debes et al. in prep), and the mineralogy and c hemistry in the 
disk wall are also d ifferent from the outer disk dThi et alJl201 it 
Mulders et al.l201 lb . A model where the region < 20 AU is dom- 
inated by small grains (> 50%), and the outer disk by large 
grains (> 99.9%) provides a good fit to both scattered light im- 
ages and SED (Fig. [8}. 

It is interesting to see that this model agrees q ualita tively 
well with the mo del presented in iBenistv et al.l (1201 Oh and 
Tatul li et alj d201 ll) : a disk wall dominated by small grains that 
produce the silicate feature, and larger grains further out. This re- 
sult appea rs not to be in agreement with the polarimetry results 
from Ouanz et al who find indeed a change in dust grain 

properties around 50... 100 AU, but with small grains further out 
and bigger grains inwards. However, it should be noted that this 
conclusion is mainly based on brightness asymmetry, which we 
have shown might not be unambiguous in deriving grain sizes. 




0.1 1.0 10.0 100.0 

X [ytim] 



Fig. 9. Observed SED of HD 100546 plus photosphere (grey 
lines). Overplotted are a series of models with selfconsistent dust 
settling with a turbulent mixing strength of a = 0.01, for grains 
with surface-to-mass-ratios of cr/m= 10 5 cm 2 /g (solid line) and 
4 • 10 3 cm 2 /g (dotted line), corresponding to compact grains of 
0.1 and 2.5 yum, respectively. 



4.2. Dust settling 

Apart from t he pre sence of m i llimet er-sized grains 

dBouwman et all 120031; IBenistv et all 1201 Ol) - which are 
most likely located near the disk midplane - the SED of HD 
100546 with its strong mid and far-infrared excess does not 
show strong signs of dust settling. A population of small grains 
in a hydrostatic disk ca n explain the SED (Do minik et al.l 
120031; iMulders et al.ll201ll) . showing that dust and gas must be 
well-mixed. In this section we explore if the large grains we find 
can be present at the disk surface. 

To find the maximum grain size that can be present at the 
disk surface, we run a s eries of models that use sel f-consistent 
settling as described in Mul ders & Dominikl d2012l) . The dust- 
gas coupling in this model depends on the dust-to-gas ratio (dtg), 
turbulent mixing strength ff tur b and dust grain surface-to-mass ra- 
tio (cr/m). For plausible values of these parameters (dtg = 0.01, 
<*turb = 0.01, we find that in a disk that has 0.005M o of gas a 
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minimum surface-to-mass ratio of 10 5 cm 2 /g is required to keep 
particles present at the disk surface (Fig . [9j . This corresponds to 
a 0.1 fim compact particle. 

To keep particles larger than 2.5 fim present at the disk sur- 
face would require a stronger dust to gas coupling. We can think 
of three ways to achieve this: 

1. A turbulent mixing strength higher than a tur \, = 0.2. 
Although such a hi gh a could be present in the upper lay- 
ers of the disk (e.g. lFromang & Nel son 2009), it would make 
HD 100546 a special case among Herbi g stars which have on 
avera ge a much lower mixing strength ( Mulders & Dominik 
120 12h . and it is beyound the scope of this paper to explore 
this further. 

2. A gas mass higher than O.1M . Although such a mass 
is not unreasonable f or a 2.4 M Herbig star (e.g. 
IWilliams & Cieza 20llT) . observations of gas emission lines 
point towards a lower rather than a higher gas mass, with a 
total mass in th e range 0.0005... 0.0 1M Q (iPanic et al]l2010t 
iThi et al]l20Tlt) . 

3. Porous grains. A more porous grain has a higher surface to 
mass ratio than a compact particle, resulting in a stronger 
dust-gast coupling. Keeping the dust mass fixed, the surface- 
to-mass ratio of an aggregate scales with filling factor ff as 



m 



4na 2 



ff- 2 / 3 . 



(4) 



because the radius scales as a — fT 1/3 . For an increase of a 
factor of >25 in the surface-to-mass ratio this gives a filling 
factor of ff< 0.01, a very fluffy aggregate. 

Aggregates can explain the presence of large particles at 
the disk surface, and p otentially the o bserved strength of the 
silicate feature as well dMin et ai]l2006l). In addition, they are 
predicted by dust g rowth experiments (Wurm & Blum 1998; 
Dom inik et al.1 2007 ) and simulations (e.g. lOrmel et alJ 12007: 
Zso m & Dullem ond 2008[). 



4.3. Phase function at intermediate scattering angles 

However, the scattering properties of complex aggregates are 
less well known than those of compact spheres. The general scat- 
tering behaviour is expected to be similar for different particle 
structure - i.e. a strong forward scattering peak due to diffrac- 
tion and a red color of the effective albedo - while the detailed 
behavior of light scattering at intermediate angles is less well 
established: the phase function for reflection/refraction is less 
well known and can be forward or backward s cattering depend - 
ing on the shape and structure of the particles dMin et al. 2010). 
Since the phase function can only be derived from observations 
at intermediate angles, fitting the brightness asymmetry requires 
including the aggregate structure. 

Full Discrete Dipole Approximation calculations are re- 
quired to calculate agregate scattering properties, which makes 
them computationally impractical to use directly in radiative 
transfer codes. The beginnings of a computationally less expen- 
sive theory for computating a ggregate opaci ties based on effec- 
tive medium theory are there dMin et alJl2006h . However, effec- 
tive medium theory has the disadvantage that the main assump- 
tion is that the constituents mixed are much smaller than the 
wavelength of incident radiation. This implies that the scattering 
efficiency of the constituents in this approximation is zero. We 
believe that in reality the constituents of the aggregates are on the 
order of a micron, which makes the scattering properties of the 



constituents an important aspect for computing the phase func- 
tion of the aggregate at intermediate scattering angles. Therfore 
their scattering properties need to be tested before they can be 
applied directly to scattered light images in the way we have 
done in this paper. 

5. Conclusion 

We have studied the effects of grain size on the colors, surface 
brightness and brightness asymmetry of scattered light images 
of circumstellar disks. We have used a 2D radiative transfer 
code that includes anisotropic scattering to model the disk of 
HD 100546 over a broad wavelength range (0.4 to 2.2 fim). Our 
conclusions are: 

- The low observed albedos of circumstellar dust can be ex- 
plained by extreme forward scattering by grains larger than 
the observing wavelength. This reduces the effective albedo 
below the practical lower limit of 0.5 for the real albedo, and 
produces gray to red disk colors. Small grains also have low 
albedos, but with very blue colors. 

- The brightness asymmetry between the front and back side 
of a circumstellar disk is not a unique indicator of grain size. 
For large grains, the forward scattering peak in the phase 
function is not observed, while the phase function in the 
observed range of angles can appear more isotropic. This 
mimics the brightness asymmetry of small grains, but with 
a lower albedo. 

- The red colors, low albedo and small brightness asymme- 
try of the protoplanetary disk around HD 100546 can be ex- 
plained by a grain size larger than 2.5... 10 micron, depend- 
ing on composition and particle structure. The presence of 
such large particles at the outer disk surface indicates they 
must be present in the form of porous aggregates, rather than 
compact particles. 

While we believe that large aggregates, rather than compact 
particles, are responsible for the observed emission, further re- 
search in characterizing the phase function of extremely forward 
scattering porous aggregates is necessary to explain their prop- 
erties on basis of scattered light images. 
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